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Dissipative preparation of a pure steady state usually involves a commutative action of a coherent and a 
dissipative dynamics on the target state. Namely, the target pure state is an eigenstate of both the coherent and 
dissipative parts of the dynamics. We show that working in the Zeno regime, i.e. for infinitely large dissipative 
coupling, one can generate a pure state by a non commutative action, in the above sense, of the coherent and 
dissipative dynamics. A corresponding Zeno regime pureness criterion is derived. We illustrate the approach, 
looking at both its theoretical and applicative aspects, in the example case of an open XXZ spin-1/2 chain, 
driven out of equilibrium by boundary reservoirs targeting different spin orientations. Using our criterion, we 
find two families of pure nonequilibrium steady states, in the Zeno regime, and calculate the dissipative strengths 
effectively needed to generate steady states which are almost indistinguishable from the target pure states. 

PACS numbers: 75.10.Pq, 03.65.Yz, 05.60.Gg, 05.70.Ln 


I. INTRODUCTION 

One indispensable pre-requisite for quantum information 
processing is preparing a given quantum state and maintain¬ 
ing it for a sufficiently long time. A promising perspective 
in generating quantum states with desired properties is of¬ 
fered by using a controlled dissipation. Instead of produc¬ 
ing a detrimental decoherent effect on the quantum system, 
the controlled dissipation can help to create and preserve the 
coherence. With the help of the controlled dissipation, one 
can prepare and maintain entangled qubit states urn, per¬ 
form universal quantum computational operations IIIH3, gen¬ 
erate and replicate entanglement between macroscopic sys¬ 
tems GMl, and store and protect quantum memory m. 
Dissipative state engineering methods are robust, since, due to 
the dissipative nature of the process, the system is driven to¬ 
wards its nonequilibrium steady state (NESS) independently 
of the initial state and of the presence of perturbations. 

Dissipative pure state engineering typically requires com¬ 
mutative actions on the target state by both the coherent and 
dissipative parts of the effective dynamics 0161EHS]. In 
other words, the target state is required to be an eigenstate 
of the Hamiltonian and of all quantum jump operators, see 
Eq. QIH. On the other hand, generic non-commutative 
coherent and dissipative actions result in a mixed steady 
state 1201 . 

In this paper, we demonstrate that by applying sufficiently 
strong dissipative couplings, one can generate steady states, 
which are arbitrarily close to pure states, for non-commutative 
dissipative and coherent dynamics. Namely, while the target 
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pure state is still required to be an eigenstate with respect to 
the quantum jump operators, it is not generically an eigen¬ 
state of the Hamiltonian. This is not in conflict with previous 
results, since the exact pure NESS is attained only in the Zeno 
limit, i.e. in the limit of infinitely strong dissipative action, 
where the NESS pureness criteria nsiiigi are not valid. 

The Zeno regime belongs nowadays to a standard toolbox 
of dissipative protocols ED- It is usually associated with an 
effect of freezing the whole quantum system, or freezing some 
degrees of freedom and accelerating some others (static Zeno 
effect, dynamic Zeno effect, anti-Zeno effect) Il22lj2^ . In the 
following, we derive a general criterion of steady-state pure¬ 
ness which applies exactly in the Zeno regime but can be used 
to generate an almost pure NESS for sufficiently strong dis¬ 
sipative couplings. We demonstrate the applicability of our 
criterion by obtaining two classes of pure stationary states 
in nonequilibrium boundary-driven Heisenberg XXZ spin 
chains, both in the critical and noncritical phases. Moreover, 
we show that, in practice, reaching the Zeno regime is not 
necessary, since applying a dissipation above a finite strength 
is sufficient to obtain pure steady states with arbitrary pre-set 
pureness. 

II. ZENO REGIME PURE NESS CRITERION 

We consider an open quantum system in contact with an 
external environment. The effective time evolution of the re¬ 
duced density matrix p of the system is described by a quan¬ 
tum master equation in the Lindblad form Il25lj27l . 

^^=-,[H,p]+TV[p], (1) 

where H is the Hamiltonian representing the coherent part 
of the evolution, E measures the strength of the dissipative 
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coupling, and I?[/5] is the Lindblad dissipator, 

V[p] = ^ U^pLi - ^ (LiL^p + pLiL^)) , (2) 

a ^ ' 

defined in terms of a set of Lindblad, or quantum jump, oper¬ 
ators, {La}- We set = 1 and J = 1, where J is a global 
energy factor which multiplies H, measuring energy in units 
of J, time in units of H/J, and L in units of J/H. A NESS is 
a fixed point solution of the dynamical Lindblad equation Q. 
We shall assume that the NESS is unique. It is easy to see that 
the NESS is a pure state, namely, pNESs(r) = |'P)('I'|, if |'I') 
is an eigenstate of the Hamiltonian and a dark state (i.e., an 
eigenstate with zero eigenvalue) with respect to all Lindblad 
operators, 

£f|«')=A|^') and = 0, for all a. (3) 

Most theoretical studies and experimental protocols rely on 
this sufficient condition ([^ for dissipatively preparing pure 
states. It often happens, however, that for the given set 
iL, {La}, no pure state satisfying the conditions ([^ lfT9l can 
be found. In those cases, it is worth formulating a less de¬ 
manding criterion by requiring pness (r) to become pure only 
in the Zeno limit E —oo. We then assume that for sufficiently 
large E, the following expansion in powers of (1/E)^ exists: 

OO 

PNESs(r) = |4')('I'|+^E-V"\ (4) 

k^l 

where the first term of the expansion p^^'> = l^k) ('I'| is a pure 
state. Inserting the time-independent state Q into Eq. Q and 
comparing the terms at different orders of E, we obtain 

I?[|4-)(^|] =0 (5) 

and the recurrence relations 

fc = 0,l,2,..., (6) 

which have the formal solution 

p("'+i) =2?-i[*[if,p(")]], fc = 0,l,2,.... (7) 

The existence of is granted if and only if 

[H, lies entirely in the subspace orthogonal to the kernel 
of V, i.e. 

PkeM[H,p<-’^^]) = 0, k = 0,1,2..., (8) 

where Pq denotes the orthogonal projector on H. In particular, 
the zeroth-order condition reads 

Pkcri5([Lf,|4')(«'|]) = 0. (9) 

Conditions Q and (|^, which, for brevity, will be named the 
Zeno regime pure NESS criterion, substitute the criterion 0 
in the limit E —oo. As we will demonstrate in the follow¬ 
ing, the Zeno regime pure NESS criterion is less restrictive 
than criterion ([^ IT^ . Moreover, satisfying Eq. Q and just 


the zeroth order necessary condition (j^ can be enough to find 
a pure NESS in the Zeno limit. By continuity, for sufficiently 
large dissipative coupling E, the actual NESS will be arbi¬ 
trarily close to the pure state. The target state I'k) is not an 
eigenstate of the Hamiltonian H ; otherwise the condition (j9]) 
becomes trivial. On the other hand, the condition Q implies 
MM that the target state is an eigenstate of the quantum 
jump operators {La}. Thus, the actions of the coherent and 
dissipative parts of the dynamics on |'I') are non-commutative, 
HLaY^) 7 ^ LaH\'^), which implies that the target pure state 
cannot be exactly prepared for any finite E. 


III. HEISENBERG SPIN CHAINS 

To test the Zeno regime pure NESS criterion, we consider 
an open XXZ Heisenberg spin chain with Hamiltonian 

^ ^ E ^ - ^)) ’ (10) 

i=i 

where A is the dimensionless anisotropy parameter measur¬ 
ing the ratio between the couplings of the Z and XY spin 
components, and a dissipator with just two Lindblad opera¬ 
tors, Li = Cl and L 2 = Cr acting locally on the “left” and 
“right” boundary spins only. The operators Cl and Cr fa¬ 
vor an alignment of the boundary spins at fc = 1 and k = N 
along the vectors Il,Ir defined by longitudinal and azimuthal 
coordinates as 

II = (sin6*L cos sin6*L sin cos 6 *l), 

Ir = {sin Or cos pR, sin Or sin pR, cos Or ). 

If IL 7 ^ Ir, then there is a boundary gradient leading to a 
NESS with nonzero current. Eor specific boundary gradients, 
the NESS of this model has been calculated analytically at 
arbitrary dissipation strength 112814311 . The explicit form of 
Cl,Cr is given in Appendix [A] where we also detail the con¬ 
tent of Eq. (j^ and the calculation of the super-operator inverse 
I)~^. In the Zeno limit, the boundary spins 1, iV are projected 
into the states described by the one-site density matrices 


PL = ]^{l PIl - CTl) , 

(11) 

Pr = ■ 

(12) 

These are single qubit pure states, tr p\ = tr plj = 
We look for a zeroth-order pure NESS, pness 
in the factorized form 

: 1. 

= \^){H 

1^') = IV'i) ® IV'2) O ■ ■ ■ \^n), 

(13) 


with IV'fc) satisfying the generalized divergence condition 

h {{lYk) ® \i’k+i)) = ftfclV’fe) ® \i’k+i) 

+ Wk) ® l^/’fc+i) - IV'fe) ® \Uk), (14) 
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where /i is a local density of the Hamiltonian ( [Io| and \Uk) 
is a local unknown vector. Substituting expression into 
Eq. (|^, we find that this is satisfied if and only if 


N-l 

“ N) = 0 . 


(15a) 

k^l 

Rk — Rk-> ^ — 1, . . 

.,A, 

(15b) 

ti Ri = tr = 0? 


(15c) 


where Rk = |C/fe)(V'fc| 
|'0fc)(C^fc-l|- 

Proof. Denoting pk 
rewrite the commutator 


|V’fe)(17fe| and Rk = \Uk-i){'fk\ - 


= |V'fc)(V'fc| and using Eq. ([^, we 
[Ef, Pness] as 


N-l 

[H, pness ] = ^ Pi 0 ■ • • 0 {Rk — Rk) ® ® Pn 

k=2 

+ 0 P2 0 • • • 0 PAf - Pi 0 • • • 0 PW-1 0 Rn 

N-l 

+ {pk - Pk) PRESS' ( 16 ) 

fc=l 

Requiring that Eq. (j^ is satisfied and taking into account that 
tr(pi) = tr(pAr) = 1 and tr(A B) = tr(A)tr(f3), we 
obtain 

The criterion ([^ would not, in the present example, provide 
any nontrivial solution; the NESS is not pure for any finite E 
and for any boundary polarization gradient. The only solution 
of Eqs. ([^ is obtained for identical boundary conditions, II = 
Ir, and fixed anisotropy, A = 1, and it corresponds to a trivial 
ferromagnetic state p = (pl)^^■ Conversely, using the Zeno 
regime pure NESS criterion, we readily find the following two 
nontrivial families of solutions. 


to meeting the boundary conditions ^pl =Q and = ^- We 
conclude that, for any twisting angle $, we have a factorized 
state which satisfies the Zeno regime conditions Q and (|^ 
only when the anisotropy assumes the values 

A(<i), m) = cos (($ + 27rm)/(A — 1)), (18) 

with m = 0,1,..., A — 2. This solution represents an 
equidistant twisting of the polarization vector in the XY plane 
along the chain with winding number m; see Eig. [T] for an 
illustration. Eor a fixed twisting angle $ and in the limit 
N —>■ oo, the set of the anisotropies ( [T8] l becomes dense in 
the interval [—1,1]. While the pure states that we obtain, are, 
by construction, dark states of the quantum jump operators, 
= 0, a = 1, 2, they are not eigenstates of the Hamilto¬ 
nian, H\'^) f AjT'). 



A. Boundary twisting in the XY plane 


Let us choose the boundary polarizations in the XY plane. 
Due to isotropy, this choice can be parametrized by a single 
angle T* between the left and right boundary polarizations, 
i.e., we can put = (1,0,0), = (cos <&, sin $, 0). Var¬ 

ious properties of the XXZ model with boundary twisting in 
the XY plane for strong and weak driving have been inves¬ 
tigated for $ = 7 r /2 and arbitrary A in ll^ 1^ . while for 
the isotropic case A = 1, the full analytic NESS for arbitrary 
r, $ has been obtained in I 


We look for a solution of Eq. (14 1 taking j'i/jfe) in the form 


IV'/c) = 


1 ( 
V 


(17) 


which corresponds to a local spin polarization Ik = 
{cos (fk, sin if k,0). As detailed in Appendix [b] such a solu¬ 
tion exists and, via Eq. (dni, is also a solution of Eq. (j^, pro¬ 
vided that (fik+i — Pk=l A = cos( 7 ). The constant 7 is 
fixed by requiring that Eq. Q is also satisfied, which amounts 


Figure 1: (Color online) Orientation, for a pure NESS in the Zeno 
regime, of the lattice spins of a driven XXZ chain arranged along 
the Z axis with boundary twisting in the XY plane for winding num¬ 
ber m = 0,1,2 (left to right). The Z labels indicate the lattice points 
along the chain. Parameters: A=12,‘l> = 7r/3. 


Since Eqs. Q and (j^ satisfied by our XY -twisting solution 
are just necessary (but not sufficient) conditions for the NESS 
in the Zeno limit to be pure, one needs an independent check 
of the pureness of the found solution. A straightforward ana¬ 
lytic computation for small system sizes N < 7 reveals that 
indeed the found NESS in the Zeno limit becomes pure ex¬ 
actly for the anisotropies ( fTS] ), with two exceptions, namely, 
$ = 0 and A = 0; see Appendix[C| Moreover, we find that no 
other pure states in the Zeno limit exist. Thus, all solutions of 
Eqs. (Bland ([^ for real-valued A are given by the factorized 
states p 3 } with anisotropy ( |T8] l. 

In Eig. 1^ we show the von Neumann entropy S — 
— tr(pNESS log2 Pness) versus the anisotropy A, in the Zeno 
limit and for finite E, obtained numerically for a system of 
four sites. In the Zeno limit, the NESS becomes pure, i.e., 
5 = 0, only at the points predicted by Eq. ( [T8] l. Eor finite 
E, the NESS is always mixed. However, at the points (I 81 , 
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and for F finite but larger and larger, pness approaches the 
respective pure states arbitrarily closely. 



Figure 2: (Color online) The von Neumann entropy S = 
— tr(plog 2 p) for p = Pnbss versus the Z-axis anisotropy A, in 
a driven XXZ chain with boundary twisting in the XY plane. The 
dotted, dot-dashed and dashed lines correspond to finite dissipative 
couplings, r = 2,5,15, respectively, while the thick line is the Zeno 
limit. Parameters: A = 4, = tt/S. 


To find out if the pure states found in the Zeno regime 
are experimentally accessible, we have numerically calcu¬ 
lated the minimal dissipation strength T^{m, N, $) required 
to reach a pure NESS within a given tolerance e, and the re¬ 
laxation time needed to establish the NESS, namely, the in¬ 
verse gap A(re)“^ of the spectrum of the Liouvillian £[•] = 
—i [H, ■] + ri2[-] at r = Fg. In practice, we define Fg as 
the dissipation strength at which the von Neumann entropy of 
the corresponding NESS becomes equal to e; see Appendix 
[D] for details. Most remarkably, we find, on the base of a 
study of small-size systems (N < 9), that the optimal (mini¬ 
mized among all the winding numbers m M) Tg decreases 
with N, making the effective “Zeno regime” more and more 
accessible as the system size increases; see Eig.|^ This some¬ 
what counter-intuitive property follows from the fact that for 
longer chains, it becomes easier to freeze the boundary spins, 
i.e. to suppress their fluctuations, so that the effective Zeno 
regime is reached earlier. In compensation, the corresponding 
relaxation time increases with N-, see Eig. However, this 
increase is only polynomial, in accordance with the general 
observation made in E). 


N 

2 3 4 5 7 10 



Figure 3: (Color online) Minimal dissipation strength Fg needed to 
reach a pure state with accuracy e = 10“^ versus system size N 
(dots, left-bottom axes). The solid line is the fit Fg = 94. 

Gap Aug of the Liouvillian spectrum at F = Fg versus system size N 
(squares, right-top axes). The solid line is the fit Apg = 9.42^“^ ®®. 
Parameters: $ = tt/3, A = cos('I>/(A’ — 1)). 


which corresponds to a local spin polarization Ik — 
(sin dk,0, cos 6k), and then restrict the found solution to meet 
Eqs. 0 and (|^. As detailed in Appendix]^ we have a NESS 
which can be of two kinds, corresponding to the orbital an¬ 
gles dk monotonically decreasing or increasing in the interval 
]0,7r[; see Eig. We never have a pure NESS with 6k = 0 ,tt 
unless in the thermodynamic limit N ^ oo. Eor finite-size 
systems and given boundary polarizations in the XZ plane, 
i.e., dk, 6ji, the NESS becomes pure in the Zeno limit and 
only for the anisotropy value, 

M0l,6r) [tan(0fl/2)/tan(0i/2)]'/('^-') 

-F i [tan(6»L/2)/ ta,n{6R/2)]^^^^~^^ . (20) 


In Eig. 1^ we show the dependence of the von Neumann en¬ 
tropy S = — tr(pNESS log2 Pness). obtained by numerically 
evaluating the NESS in a system of four sites for different val¬ 
ues of the anisotropy. In the Zeno limit, the NESS becomes 
pure, i.e.. S' = 0, at the point predicted by Eq. (20 1 . 


B. Boundary twisting in the XZ plane 

Next we orient the boundary polarization in the XZ plane, 
Ir = {sin6k,0, cos 6l), Ir = {sin6R,0,cos6R). As before, 
we first solve Eq. ([l4]l, now taking l-ijjk) in the form 


IV'/c) = 


f cos % 

isin^ 


(19) 


IV. CONCLUSIONS 

To summarize, we have formulated a criterion for a 
nonequilibrium steady state of an open quantum system to be 
pure, in the Zeno limit, i.e., for asymptotically large dissipa¬ 
tive coupling. The criterion is specified by Eqs. (|^ and ([^. 
Zeno-limit pure states are not reachable, in a strict mathemat¬ 
ical sense, for any finite dissipative coupling. However, by 
applying a finite but large enough dissipative coupling, one 
can generate pure NESSs with arbitrary precision. 
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Figure 4: (Color online) As in Fig. for boundary twisting in the 
X Z plane and in correspondence to solutions with decreasing (left) 
and increasing (right) orbital angles. Parameters: N = 7, 6l ~ it/2, 
Or = tt/IO (left), 6r = TT — tt/IO (right). 



Figure 5: (Color online) The von Neumann entropy S = 

— tr(plog 2 p) for p = Pness versus the Z-axis anisotropy A, in a 
driven XXZ chain with boundary twisting in the XZ plane. Param¬ 
eters: N = A, 6l ~ TT12, Or = arctan(l/(15x/3))- The NESS is 
pure for a single value of the anisotropy, A = 2, given by Eq. po^ . 
The dot-dashed and dashed lines are obtained for finite dissipative 
couplings, r = 30,370, respectively, while the thick line is the Zeno 
limit. 


Using our criterion, in the Zeno regime we find two fami¬ 
lies of pure NESSs for the driven quantum XX Z spin chain 
with boundary twisting in the XY or XZ plane, for values 
of the Z-axis anisotropy given by Eqs. ( [TS] ) and ( [20| i, respec¬ 
tively. The criterion can be straightforwardly applied to gen¬ 
erate pure steady states in other nonequilibrium quantum sys¬ 
tems. 

Our approach opens an interesting perspective in dissipa¬ 
tive engineering of pure states. If, for given resources, prepar¬ 


ing an exact pure steady state at finite dissipative strength is 
impossible (as it happens in our example of driven XXZ 
chain), it may still be possible to generate a pure state in the 
Zeno limit. In practice, this means that at finite dissipation, 
a slightly mixed state will be produced, which, however, be¬ 
comes infinitesimally close to a pure state as the dissipation is 
increased. The effective coupling needed to reach the “Zeno 
regime” depends on the chosen measure of pureness and the 
required precision and must be estimated in each case sepa¬ 
rately. In the example of the driven XX Z model considered 
here, the effective Zeno regime is reached at very moderate 
dissipative couplings. 
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Appendix A: Inverse of the Lindblad dissipator and secular 
conditions 

The Lindblad operators Li = Cl, L 2 = Cr have the form 

Li = [(cos^L cos<pl)o'i + (cosOl sm(pL)o'f 

- (sin + icrf(-smipL) + iaf(cos(pL)]/2, 

L 2 = [{cos 9n cos (fR)a^ + (cos sin 

- {sin 9 R)a^^ + ia%{-sin If r) + {cos fR)]/2. 

The dissipator, V = Vl + 'Dr, is the sum of the left and right 
dissipators 

Dl[] = l,-l\-^{l\l,,-}, 

Dr[-]=L2-LI-^{LIL2,-}, 

which are linear super-operators acting locally on a sin¬ 
gle qubit. The eigenbasis of the eigenproblem 

'^rI^r] = ^a/)R is 

= {2pR, 2pR -I,- sinipRa"^ + cos fRa^, 

cos9r{cos fRcr^ +sinipRa'‘^) — sinfl^cr^} 

with the respective eigenvalues 

{Aq} = {0, -1,--,--}. 

Here I is a 2x2 unit matrix, are the Pauli matrices, 

and pr is the targeted spin orientation at the right boundary. 
Analogously, the eigenbasis and eigenvalues of the eigenprob¬ 
lem are 

= {2pL,2pL - I, - sin ipLcr'' + cos ipLcr^, 
cos 9 L {cos ifLcr^ + sinipLcr^) — sin 
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where pi is the targeted spin orientation at the left boundary. 
Since the bases {(!)%} and {(j)^} are complete, any matrix x 
acting in the appropriate Hilbert space can be expanded as 


where h is the local density of the Hxxz Hamiltonian, 

/i = ^ (g, + a" (g) a" + A (g) - /)), 

and \Uk) is a local unknown vector, 


4 4 

X = (j)R- 

a—1 /3—1 


(Al) 


Indeed, let us introduce two complementary bases 

{V’L.fl} = {I/‘^,PL,R - I, i-sin(pL,Rcr^ + COS ipL,RCr'^)/2, 
(cos 9l,r{cos (PL,R 0 -'' + sin tpL.Rcr'^) 

- sin6>i_j^cr^)/2}, 

trace-orthonormal to the {(p^ ^j}, namely, tr('!/)J^^) = Sa^y 

and = 5py. Then, the coefficients of the expansion 

are given by 

Xpa = tri,Ar((^^ g g 


where tri at denotes the trace taken with respect to the first- 
and the last-sp in sp aces only. On the other hand, in terms of 
the expansion (All, the super-operator inverse 'D~^ = (T’l + 
'Dr)~^ is simply 


{Vr+Vr)-\x]=Y1 


Oc,^ 


Act + P/3 


< g X/3ct g (pR- 


\Uk) = 


Uk 

Vk 


Equation (B1 1 is an overdetermined system of equations for 
IJ,k,Uk,Vk- The system does not admit a solution unless the 
Z-anisotropy parameter takes the value A = cos((pfe+i — (pk), 
which is possible only if the difference between any two con¬ 
secutive angles along the chain is kept constant, (pk+i — Pk = 
7 . In this case, we have 


Pk+l 


Pfc = 4 sin 
Uk = —iv/ 2 sin 




Pk 


4 

Pk+l 


Pk Pk+l - Pk _i:fk±i±fk 
cos „ — —e 4 


Pk+l 


Pk Pk+l - Pk ^ Vk + l + Vk 
-^cos^-^- —e 4 

4 2 


Vk = iv/^sin 
From the above solution, we compute 


Rk = \Uk){ipk\ - \'Pk){Uk\ = 1^0^ sin((pfc+i - Pk), 

~ % 

Rk = \Uk-i){i>k\ - \'ipk){Uk-i\ = -cr^ sin((pfc - Pk-i). 


The above sum contains a singular term with a = /? = 1, 
because Ai -f pi = 0. To eliminate this singularity, one must 
require Xn = tri Ar X = 0^ which is equivalent to the secular 
condition 

■Pkcr-D(x) = 0, 

where Pq denotes the orthogonal projector on H. 

We conclude that the existence of 
at order fc = 0 , 1 , 2 ,... is granted if and only if 

tri,Ar([Tr,p«]) =0. (A2) 

Appendix B: Stationary states with boundary twisting in the 
XY plane 


It is straightforward to check that the system of equations. 


N-1 

i^^k-p*k) = 0 , 

(B2a) 

k=l 

Rk = Rk, fc = 1 ,..., A, 

(B2b) 

tri?i = tr i?Ar = 0 , 

(B2c) 


is then satisfied. Since Eq. ( |B2| i has been demonstrated to 
be equivalent to PkGrX>([^^, l^)('I'|]) = 0, we conclude that 
the found solution meets the necessary condition of our Zeno 
regime pure NESS criterion. To meet the other condition, 
namely, I?[|T')(T'|] = 0 , we just need to satisfy the bound¬ 
ary conditions pi = 0 and (pjv = $. This is accomplished by 
choosing 


Assuming p(°) = |'I')(T'| in the factorized form 

| 4 ') = IV'i) g IV'2) g • • • |' 0 jv), 


with 


I'fpk) 


V2 V 


/c = 1, 




7 = 7 ( 4 ), m) = (4> -f 2'Km)/{N — 1), 
with TO = 0 , 1 ,..., A — 2 . 


Appendix C: Analytic calculation of the Zeno NESS for small 
sizes and boundary twisting in the XY plane 


we look for a solution of the generalized divergence condition 

h {\'ipk) g IV'fe+i)) = Pk\i’k) g IV’fe-i-i) 

+ \Uk) g \tpk+i) - \tpk) g \Uk), 

(Bl) 


Solving the secular conditions (j^ for fc = 0,1, we compute 
the analytic form of pness{N, A, $) for small system sizes 
N, and calculate the pureness parameter, 

/(A, A, 4)) = trlpxESsiN, A, 4>)^] - 1. 
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Note that the state pness is pure if and only if / = 0. We 
obtain 


/(3,A,$) 

/(4,A,$) 

/(5,A,$) 


(-r 2 (A) + cos($))^ 
2(2A2 + cos($) + l)^’ 

(-T 3 (A)+cos($)f g, 
(-T 4 (A)+cos($)f 


where T„(cosx) = cos(nx) are Chebyshev polynomials of 
the first kind, 

T2ix) = -1 + x^, 

Taix) = —3x + ix^, 

T^{x) = 1 — 8x^ + Scc^, 


Zeno NESS is pure only for A = 1 for N even, which cor¬ 
responds to TO = 0 in Eq. ( pl , and A = ±1 for N odd, 
corresponding to to = 0, {N — l)/2 in Eq. (Cl i. The A = 1 
solution represents a trivial factorized state with all spins po¬ 
larized in the X direction. This homogeneous state remains a 
NESS for any finite value of E. 

Exception (b) A = 0. Whenever among the critical 
anisotropy values (Cli, a free fermion point A = 0 ap¬ 
pears, the respective NESS at A = 0 is not a pure state, 
but a fully mixed state, apart from the boundaries, pness = 
Pl®{I The peculiarity of this exception results 
from the fact that the Zeno limit E —^ oo and the free-fermion 
limit A —?> 0 do not commute, with the reason being the exis¬ 
tence of an extra symmetry of the NESS at A —?> 0; see 
for an elaborate treatment of a $ = 7r/2 case. 


Appendix D: Minimal dissipation strength 


and, for instance, 

P 4 =16A6 -f 56A4 -f 4 (lOA^ -f 3) A cos($) 
-I- 30A^ -I- 3 cos(2$) -I- 3, 


Qi =2 (I 6 A® -f 12A^ + 4 (5A2 -f 2) a cos(4>) 

-f 14A2-pcos(2$)-f l)^ 

Substituting A = cos 7 into the above expressions for /, 
we obtain 


/(3,cos7, 4>) 


( 003 ( 27 ) — cos(<i)))^ 
2(cos(27) -f cos($) -I- 2)2 ’ 


Eor system sizes 3 < iV < 9, we have numerically 
calculated pness( r)^ namely the NESS of the Liouvillian 
—i[H, •] + EI?[-], where H is the Hamiltonian of the XXZ 
model and V = 1)^ -f 27/j is the dissipator described in Ap¬ 
pendix]^ for several hnite dissipation strengths E. In the case 
of boundary twisting in the XY plane, the von Neumann en¬ 
tropy S'(E) = -tr(pNESs(r)log2PNESs(r)) Corresponding 
to the NESS obtained for to = 0 and $ = 7 r /3 is plotted in 
Fig- § as a function of E. We see that, for any N, S'(E) de¬ 
creases monotonously to 0 by increasing E, approximately as 
E“^ for E large. As is natural, we define the minimal dissi¬ 
pation strength Ee(TO, N, $) required to reach a pure NESS 
within a given tolerance e as the unique solution of 

-S'(Ee) = - tr [pNESsi^e) log 2 PNESsi^e)] = C, 

which is plotted in Eig. 3 of the main text. 


/(4,cos7, 4>) = - (cos(37) - cos(4>))^ 
/(5,cos7, $) = - (cos(47) - cos(4>))^ 


Extrapolating for arbitrary A^, we get cos((7V— 1 ) 7 ) = cos(<i)) 
as a pure NESS condition, yielding 

7($,to) = (4) -f 2TTm)/{N — 1), to = 0,1, ...,N — 2, 

as well as 

A($,to) = cos(7($,to)) = cos ( I- (Cl) 



Inde pend ently, we verify that the anisotropy values given by 
Eq. ( |C1| exhaust the solutions of the equations /(A, A, T*) = 
0 for A being real. Points of non-analyticity of the functions 
f{N, A, 4)), e.g. /(3, 0, tt) = 0/0, correspond to exceptions 
and need to be analyzed separately. 

Exception (a) 4) = 0. This case corresponds to a full bound¬ 
ary alignment, i.e. to the absence of a boundary gradient. The 


Figure 6: (Color online) The von Neumann entropy 5(r) = 
— tr[pNESs(r) log 2 PNESs(r)] Corresponding to the NESS calcu¬ 
lated for m = 0 and 4> = vr/S in the case of boundary twisting 
in the XY plane and for system size N from 3 (top) to 9 (bottom). 
The dashed line is the eye-guide S while the horizontal solid 

line is the tolerance value S = 10“®. 
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Appendix E: Stationary states with boundary twisting in the with 
XZ plane 


Here, we solve Eq. (14 1 taking in the form 


bk = (cos6>fc+i - cos6»fc)/sin6»fe+i, 
h = {cosdk - cos6»fc_i)/sin6»fc_i. 


IV'fc) 

It is convenient to denote 



Kk = (tan(6»fc+i/2)) / (tan(6>fc/2)). 


We find that if the Z-anisotropy has the value 2A = Kk + ^j.^, 
Eq. ([T4|) has the solution 


fifc = 


Uk = 


Vk = 


(cos - cos 3 

sin 9k sin 6k+i 

cos — cos ^ 2 + ^fc+i 

sin ^ sin 

sin - sin ^ ^ Ok + 6k+i 

sin 


cos 


cos ^ cos 


Erom the above solution, we compute 


Rk = \Uk){tpk\ - = ihcr^, 

Rk = \Uk-i){ilJk\ - \'il’k){Uk-i\ = ibkcr^, 


Conditions (151, i.e., PkcrX>([^^) |'I')('I'|]) = 0^ are thus sat¬ 
isfied if Kk + k'^^ = 2A and bk = bk for all k. The latter 
condition after some algebra gives 


-1 


-1 


+ kik-l — 


(El) 


Notice that since Kk are real numbers, | A| > 1. There are two 
independent solutions of Eq. namely, Kk = z±, where 
Z± = A. i: y/ A^ — 1 are the roots of the quadratic equation 
z -I- 1/z = A. To meet the condition I?[|T')(T'|] = 0, we 
require di = 9l and 9^ = 9a. 

The solutions with Kk = z± describe orbital angles 9k 
monotonically decreasing or increasing in the interval ]0,7r[. 
Note that we never have a pure NESS with 9k — 0, tt, un¬ 
less in the thermodynamic limit iV —^ oo. In conclusion, 
for finite-size systems and given boundary polarizations in 
the XZ plane, we have one NESS in correspondence to the 
anisotropy value. 


A{9l,9r) [tan(6'fl/2)/tan(6'L/2)]^-i 

-f ^ [tan(6»L/2)/tan(6»ii/2)]^ . (E2) 
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